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NOMENCLATURE 


a  acceleration  of  a  falling  sea  spray  droplet 

<?w  c’s<it,s/t’sat'  water  activity  of  an  aqueous  solution 

B  function  of  in  the  spray  generation  function  of  Monahan  et  a).  (1986) 

c  molar  concentration  of  an  aqueous  solution 

Cr,  drag  coefficient  of  a  spray  droplet 

cp<1  1 .006  x  103  J  kg'1  UC  the  specific  heat  of  air  at  constant  pressure 

c  4.0  x  103  J  kg'1  °C~\  the  specific  heat  of  a  spray  droplet  (i.e.,  of  sea  water)  at 

constant  pressure 

dF/dr  spray  generation  function;  number  of  droplets  produced  at  the  sea  surface 
per  unit  surface  area  per  unit  time  per  unit  increment  in  initial  droplet  ra¬ 
dius 

dF/dr^Q  spray  generation  function  when  all  droplets  are  brought  to  a  relative  humidi¬ 
ty  of  80% 

dF0/drm  contribution  to  dF/drS0  from  bursting  bubbles;  droplet  sizes  referenced  to  a 
standard  humidity  of  80% 

dFx/drm  contribution  to  dF/drS0  from  spume  droplets;  droplet  sizes  referenced  to  a 
standard  humidity  of  80% 

Dh  1 .4  x  10“7  m2/ s,  diffusivity  of  heat  in  sea  water  at  0°C 

Ds  6.8  x  10'10  m2/s,  diffusivity  of  NaCl  in  sea  water  at  0°C 

Dvv  bulk  molecular  diffusivity  of  water  vapor  in  air 

Dw'  water  vapor  diffusivity  in  air  modified  for  noncontinuum  effects 

csat  pressure  of  water  vapor  in  saturation  with  a  plane  surface  of  pure  water 

l’sat  s  pressure  of  water  vapor  in  saturation  with  the  plane  surface  of  an  aqueous 

solution 

/  RH/100,  fractional  relative  humidity 

fx  rate  of  change  of  the  variable  x  with  time 

/0  1  -  0.000537  S,  fractional  relative  humidity  at  the  sea  surface 

S  9.82  m/s,  acceleration  due  to  gravity 

gift  function  in  Fitzgerald's  (1975)  model  of  droplet  size 

ka  thermal  conductivity  of  air 

ka'  thermal  conductivity  of  air  modified  for  noncontinuum  effects 

Lv  latent  heat  of  vaporization  of  water 

m  molality  of  a  solution 

Ma  28.9644  x  10'3  kg/mol,  molecular  weight  of  air 

wd0  initial  mass  of  a  sea  spray  droplet 

Ms  58.443  x  10'3  kg/mol,  molecular  weight  of  sodium  chloride 

»is  mass  of  salt  contained  by  a  spray  droplet 

Mw  18.0160  x  10'3  kg/mol,  molecular  weight  of  water 

mw  mass  of  water  contained  by  a  spray  droplet 

vi ^  initial  mass  of  water  contained  by  a  spray  droplet 


v 


P  atmospheric  pressure 

P0  1013.25  hPa,  standard  atmospheric  pressure 

R  8.31441  J  mol'1  K'1,  universal  gas  constant 

r  radius  of  a  sea  spray  droplet 

rd  radius  of  the  (dry)  salt  particle  left  when  a  spray  droplet  evaporates 

radius  of  a  droplet  in  equilibrium  with  its  environment 
r0  initial  radius  of  a  sea  spray  droplet 

r80  radius  of  a  droplet  in  equilibrium  with  an  environment  at  80% 

relative  humidity 

Re  Reynolds  number  of  a  droplet 

RH  100  /,  relative  humidity 

S  surface  salinity  of  the  ocean 

Sd  instantaneous  salinity  of  a  spray  droplet 

Sv*  function  of  temperature  used  in  finding  the  apparent  molal  volume  of  an 

aqueous  solution 

T  instantaneous  droplet  temperature  or  any  general  temperature 

t  time 

Ta  air  temperature 

Teq  temperature  of  a  droplet  in  thermal  equilibrium  with  its  environment 

Tw  surface  temperature  of  the  ocean 

T0  273.15  K  (=  0°C),  the  ice  point 

U  wind  speed  at  a  reference  height  of  10  m 

u  instantaneous  speed  at  which  a  droplet  is  falling 

uf  terminal  fall  speed  of  a  droplet 

va  apparent  molal  volume  of  an  aqueous  solution 

vM  apparent  molal  volume  of  a  solution  at  infinite  dilution 

x  either  r  or  T 

xn  discrete  value  of  x  at  time  step  n 

if  function  that  gives  the  effects  of  curvature  and  dissolved  salt  on  the  vapor 

pressure  at  the  surface  of  a  spray  droplet 
a  function  of  /  in  Fitzgerald's  (1975)  model  of  droplet  size 

ac  0.036,  a  dimensionless  constant  used  in  computing  the  modified  water  vapor 

diffusivity 

a,-  0.7,  a  dimensionless  constant  used  in  computing  the  modified  thermal  con¬ 

ductivity 

P  function  of  /  in  Fitzgerald's  (1975)  model  of  droplet  size 

8  (T/Ta)  - 1,  see  eq  2 

AT  2.17  x  10"7  m,  a  length  scale  used  in  computing  the  modified  thermal  conduc¬ 

tivity 

Aw  8  x  10-8  m,  a  length  scale  used  in  computing  the  modified  water  vapor  diffu¬ 

sivity 

At  time  step 

X  mean  free  path  of  air  molecules 

v  2,  the  number  of  ions  into  which  a  sodium  chloride  mole  ale  dissociates 

va  kinematic  viscosity  of  air 

pa  density  of  air 

ps  solution  density  of  a  spray  droplet 

p^  initial  solution  density  of  a  spray  droplet 

pv  ambient  water  vapor  density 


vi 
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pvr  water  vapor  density  at  the  surface  of  a  spray  droplet  of  radius  r 

pw  density  of  pure  water 

Gs  surface  tension  of  a  plane,  aqueous  solution 

ow  surface  tension  of  a  plane,  pure  water  solution 

Tf  time  required  for  a  spray  droplet  with  terminal  fall  speed  //,  to  fall  1  m  in  still 

air 

Th  time  required  for  the  interior  of  a  spray  droplet  to  be  everywhere  within  4% 

of  an  impulsively  applied  surface  temperature 
Tr  time  required  for  a  spray  droplet  to  come  to  within  t>_1  of  its  equilibrium  radi¬ 

us 

tT  time  required  for  a  spray  droplet  to  come  to  within  t,_i  of  its  equilibrium  tem¬ 

perature 

ttc  tt  va^ue  estimated  from  eq  74 

ts  upper  bound  on  the  time  required  for  excess  salt  at  a  droplet  surface  to  dif¬ 

fuse  throughout  the  interior  of  the  droplet 
time  derivative  of  either  r  or  T  at  time  step  n 
(<t>n  +  <Cn+1)/2 

practical  osmotic  coefficient  of  sodium  chloride  dissolved  in  water 


Thermal  And  Size  Evolution  Of  Sea  Spray  Droplets 


EDGAR  L  ANDREAS 


INTRODUCTION 

As  the  wind  speed  at  a  height  of  10  m  reaches  3-4  m/s,  the  ocean  surface  becomes  dis¬ 
rupted — whitecapping  begins  (Monahan  1971,  Wu  1979,  Monahan  and  O  Muircheartaigh 
1981,  Monahan  et  al.  1983).  These  whitecaps  form  when  clouds  of  air  bubbles  that  were 
trapped  by  breaking  waves  rise  to  the  surface  and  burst.  The  bursting  bubbles,  in  turn, 
throw  sea  spray  droplets  into  the  air.  When  the  wind  speed  reaches  9  m/s,  a  second 
mechanism  for  generating  sea  spray  comes  into  play.  Now  the  wind  is  strong  enough  to 
tear  off  the  wave  crests  and  to  propel  spray  directly  into  the  air  (Monahan  et  al.  1986). 

The  effects  attributable  to  this  spray  are  manifold.  For  example,  spray  droplets  can 
carry  electrostatic  charges  away  from  the  sea  surface;  the  resulting  atmospheric  electrifi¬ 
cation  can  manifest  elsewhere  as  a  thunderstorm  (Blanchard  1963,  Roll  1965).  Spray  drop¬ 
lets  can  transport  organic  matter,  such  as  bacteria,  from  the  sea  to  the  air  (Blanchard 
1983).  Sea  spray  droplets  are  saline;  when  a  droplet  evaporates,  it  leaves  behind  a  micro¬ 
scopic  sea-salt  particle  that  the  wind  can  easily  carry  long  distances.  These  sea-salt  parti¬ 
cles  are  common  cloud  condensation  nuclei  (Woodcock  1952,  Mason  1957,  Cipriano  et  al. 
1987)  and  are  also  implicated  as  condensation  nuclei  for  marine  fogs  (Woodcock  1978).  If 
the  sea-salt  aerosol  encounters  a  more  humid  environment,  the  particles  can  become  re¬ 
wetted  to  form  marine  haze,  which  affects  the  extinction  of  electromagnetic  waves  propa¬ 
gating  in  the  marine  boundary  layer  (Schacher  et  al.  1981). 

As  the  wind  speed  increases,  the  number  of  spray  droplets  produced  increases  by  sev¬ 
eral  orders  of  magnitude.  As  a  result,  in  high  winds,  sea  spray  droplets — which  start  at 
the  temperature  of  the  ocean  surface — effectively  increase  the  oceanic  surface  area.  Spray 
thus  has  the  potential  for  enhancing  the  transfer  of  all  constituents  that  are  typically  ex¬ 
changed  at  the  air-sea  interface.  The  most  important  of  these  constituents  are  sensible 
heat,  latent  heat  (moisture)  and  gases.  Several  authors  have  already  looked  at  the  impact 
of  sea  spray  on  heat  and  moisture  transfer  across  the  air-sea  interface  (Bortkovskii  1973, 
1987,  Borisenkov  1974,  Wu  1974,  Ling  and  Kao  1976,  Ling  et  al.  1978,  Wang  and  Street 
1978,  Street  et  al.  1978,  Mestayer  and  Lefauconnier  1988)  but  have  reached  no  consensus 
as  to  if  or  when  the  spray  transfer  is  important  or  on  how  to  par;  meterize  its  magnitude. 
In  general,  most  agree  that  in  high  winds  (>  15  m/s),  the  sea  spray  must  enhance  the  in¬ 
terfacial  transfers  of  heat  and  moisture,  but  the  magnitude  of  the  effect  is  unknown. 

I,  too,  am  interested  in  what  impact  sea  spray  has  on  air-sea  heat  and  moisture  trans¬ 
fer,  especially  over  the  high  latitude  ocean.  Here,  enormous  sensible  and  latent  heat  flux¬ 
es  occur  when  air  flowing  off  the  cold  pack  ice  encounters  the  relatively  warm  ocean. 
These  fluxes  provide  energy  for  the  generation  and  maintenance  of  polar  lows  (Rasmus- 


sen  1985,  Shapiro  et  al.  1987),  small  synoptic  or  subsynoptic-scale  cyclones  that  form  in 
the  ice  edge  region  and  intensify  as  they  move  over  the  relatively  warm  waters  of  the 
open  ocean  (Kellogg  and  Twitchell  1986,  Rasmussen  and  Lystad  1987).  The  high  winds 
associated  with  polar  lows,  of  course,  produce  sea  spray;  tills  spray  may,  in  turn,  feed 
heat  and  moisture  back  into  the  system  and,  thus,  further  intensify  it.  My  long-range  ob¬ 
jective  is  to  quantify  the  magnitude  of  this  feedback. 

The  thermodynamics  of  sea  spray  droplets  at  high  latitudes  is  conceptually  simple 
(Fig.  1).  When  a  cold  wind  blowing  off  the  pack  ice  encounters  the  relatively  warm  ocean, 
it  fosters  intense,  upward  turbulent  fluxes  of  sensible  and  latent  heat  at  the  air-sea  inter¬ 
face.  By  disrupting  the  sea  surface,  the  wind  also  produces  spray  droplets.  Since  these 
droplets  must  initially  have  the  same  temperature  as  the  ocean  surface,  they  transfer  ad¬ 
ditional  sensible  heat  to  the  air  as  they  cool.  The  droplets  also  begin  evaporating  because 
the  cold  air  generally  has  low  absolute  humidity.  This  evaporation  must  extract  latent 
heat  from  the  air  but  also  provides  water  vapor  to  the  atmospheric  boundary  layer,  there¬ 
by  increasing  its  humidity  and  providing  vapor  that  immediately  condenses  to  fog  or  can 
easily  rise  to  higher  altitudes  to  form  clouds. 

Although  conceptually  simple,  sea  spray  thermodynamics  is  parametrically  complex. 
The  number  and  size  of  the  spray  droplets  produced  depend  on  the  wind  speed  and,  to 
some  extent,  on  the  air  and  sea  temperatures.  The  evolution  of  the  temperature  and  size 
(radius  or  moisture  content)  of  a  spray  droplet  depends  on  its  size,  on  the  air  and  sea 
temperatures,  on  the  ambient  humidity,  on  the  surface  salinity  of  the  ocean,  and  on  the 
wind  speed  and  the  associated  turbulent  wind  field  in  the  atmospheric  surface  layer.  In 
addition,  in  the  ice  edge  region,  the  environment  in  which  the  spray  droplets  evolve  is 
not  horizontally  homogeneous  or  in  steady  state.  As  the  air  flows  off  the  ice  and  over  the 
open  ocean,  oceanic  heat  and  moisture  continually  warm  and  moisten  it. 

To  simplify  the  analysis  for  this  initial  report,  I  will  focus  on  the  thermodynamics  of 
individual  sea  spray  droplets  in  a  steady-state  environment.  Pruppacher  and  Klett  (1978) 
developed  equations  for  the  rates  of  change  of  the  temperature  and  radius  of  cloud  drop¬ 
lets.  1  use,  basically,  their  equations  to  investigate  the  thermal  and  size  evolution  of  saline 
sea  spray  droplets.  Earlier,  Andreas  et  al.  (1981)  had  also  used  Pruppacher  and  Klett's 
equations  to  model  the  evolution  of  condensate  droplets.  Ultimately,  I  parameterize  the 
evolution  of  sea  spray  droplets  in  terms  of  e-folding  times — the  times  necessary  for  a 
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Figure  1.  Conceptual  model  of  heat  and  moisture 
transfer  associated  with  sea  spray. 
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droplet  to  come  to  within  e~x  of  its  equilibrium  temperature  and  equilibrium  radius  for 
the  given  ambient  conditions.  On  comparing  these  times  with  typical  Stokes  fall  rates,  we 
get  some  idea  of  the  size  range  of  spray  droplets  that  can  most  effectively  participate  in 
the  heat  transfer.  The  spray  generation  function  proposed  by  Monahan  et  al.  (1986)  pro¬ 
vides  the  necessary  information  on  the  size  of  spray  droplets  and  on  the  rate  at  which 
they  are  produced. 


SIZE  EVOLUTION  EQUATION 


In  their  Chapter  13,  Pruppacher  and  Klett  (1978)  developed  a  rate  equation  for  the 
change  in  radius  of  an  aqueous  solution  droplet.  Their  development  culminates  with  eq 
(13-26),  which  I  repeat  here 


rdr  L _ i_ 

dt  p  sRTa  f  1+8 
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where  r 

= 

instantaneous  radius  of  the  spray  droplet  at  time  t 

Ps 

= 

density  of  the  droplet 

Pw 

= 

density  of  pure  water 

ms 

= 

mass  of  NaCl  in  the  droplet 

D  ’ 

w 

= 

modified  molecular  diffusivity  of  water  vapor  in  air 

R 

= 

universal  gas  constant 

K 

= 

latent  heat  of  vaporization  of  water 

Mw 

= 

molecular  weight  of  water 

Ms 

= 

molecular  weight  of  NaCl 

T* 

= 

ambient  air  temperature 

f 

= 

fractional  relative  humidity  (e.g.,  if  RH  =  80%,  /=  0.80) 

— 

saturation  vapor  pressure  over  a  pure,  flat  water  surface  with  tempera- 
tureTa 

= 

surface  tension  of  a  flat  surface  with  the  same  salinity  and  temperature 
as  the  spray  droplet 

= 

practical  osmotic  coefficient  of  the  droplet 

V 

And 

nr 

total  number  of  ions  into  which  a  salt  molecule  in  the  droplet  dissociates. 

1  j..  /^\ 

d=JL-l=l^  rdL  (2) 

Ta  Taka'  dt 


where  T  is  the  instantaneous  droplet  temperature  (assumed  uniform),  and  ka'  is  the  modi¬ 
fied  thermal  conductivity  of  air.  Although  sea  water  contains  measurable  amounts  of 
many  dissolved  salts,  sodium  chloride,  NaCl,  is  by  far  the  main  component  (Neumann 
and  Pierson  1966).  Henceforth,  I  assume  that  sea  spray  droplets  contain  only  NaCl  as  a 
solute;  v  is  thus  2. 
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The  Nomenclature  lists  all  of  these  symbols  and  gives  the  values  of  the  ones  that  are 
constants.  As  we  go  on,  the  meanings  of  some  of  the  less  familiar  ones  should  become 
clearer. 

Before  simplifying  eq  1,  let  me  interpret  the  terms  so  we  can  see  what  affects  the  mois¬ 
ture  content  of  a  spray  droplet.  In  eq  1,  term  I,  the  ambient  humidity,  is  the  main  term 
driving  the  diffusion  of  water  vapor  toward  or  away  from  the  droplet.  This,  compared 
with  the  humidity  at  the  exact  surface  of  the  droplet,  dictates  in  which  direction  the  va¬ 
por  diffuses.  Terms  II,  III  and  IV  predict  the  humidity  at  the  droplet  surface  and  show 
that  three  distinct  phenomena  affect  it.  Term  II  results  because  the  droplet  is  not  at  the 
ambient  air  temperature;  therefore,  its  surface  vapor  pressure  is  computed  at  T  rather 
that  at  Ta.  Once  a  droplet  reaches  thermal  equilibrium,  if  the  droplet  is  evaporating,  it 
must  be  cooler  than  the  air  because  the  vapor  is  carrying  away  latent  heat.  Conversely,  if 
vapor  is  condensing  on  the  droplet,  it  must  be  warmer  than  the  ambient  air.  Term  III  re¬ 
flects  the  effects  of  surface  tension  and  droplet  curvature.  An  increase  in  this  term  in¬ 
creases  the  vapor  pressure  at  the  droplet  surface  and  thereby  enhances  evaporation  but 
retards  condensation.  Notice  the  r"1  dependence  in  term  III — the  term  becomes  increas¬ 
ingly  important  as  the  droplet  radius  decreases.  Term  IV  reflects  the  fact  that  salts  dis¬ 
solved  in  the  droplet  depress  its  surface  vapor  pressure. 

Since  1 5  I  is  rarely  larger  than  0.1,  in  eq  1  we  can  approximate 
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Thus,  eq  1  becomes 


r  dr 
dt 


Dw'  (Ta) 


Ps  RTa 


f  -  (1-5)  fl  + 
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Lv  Mw  5 
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exp  (y) 
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where 


_  2Mwgs  _  v  <&sms (Mw Ms)  (7) 

RTa  pwr  (4  7t  r3  ps  /3)  -  ms 

Since  surface  tension  and  solu'  effects  are  small  for  many  spray  droplets  in  the  size 
range  we  are  considering  (Pruppacher  and  Klett  1978,  p.  420),  y  is  small.  Therefore,  we 
also  approximate 

exp(y)  =  1  +  y.  (8) 
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Substituting  this  in  eq  6,  multiplying  out  the  products  and  discarding  terms  contain¬ 
ing  second  and  higher  order  terms  (i.e.,  82,  81/  and  S2i /),  we  obtain 


r  dr_  —  Pw  Myv^sat 
dt  psR  7a 


(/  -  1)  +5(l 

\  R  7 J 


(9) 


Because  8  contains  r(dr/dt)  (see  eq  2),  we  can  solve  eq  9  for  the  rate  of  change  of  the 
spray  droplet  radius,  dr/dt.  The  result  is 

dr  =  [(/-l)-y]r-] 

dt  PsRT*  +  f-vPs  Mw  _  (10) 

Dw'  Mwesal  (Ta)  TAka'\RTA  I 


Notice,  implicit  in  this  result  is  the  condition  for  a  droplet  to  be  in  moisture  equilibri¬ 
um  with  its  environment.  This  occurs  when 


/-!  =  •/• 


(11) 


That  is,  the  relative  humidity  is  dominant  in  determining  the  equilibrium  radius.  But  the 
mass  of  the  salt  contained  in  the  droplet  the  air  temperature  Ta,  and  the  equilibrium 
droplet  temperature,  which  affects  pw,  ps  and  os,  are  also  important. 

In  the  Numerical  Methods  section,  I  will  present  a  numerical  technique  for  solving  this 
equation  for  the  droplet  radius  as  a  function  of  time.  In  the  remainder  of  this  section,  I 
will  discuss  the  equations  I  use  to  compute  the  physical  and  chemical  variables  necessary 
to  employ  eq  10. 

The  ambient  environmental  parameters  that  I  specify  before  solving  eq  10  are  the  air 
temperature  7a,  the  relative  humidity  RH  (=  100  /),  the  sea  surface  temperature  Tw,  the 
sea  surface  salinity  S  and  the  atmospheric  pressure  P.  I  also  specify  the  initial  radius  of 
the  spray  droplet  rg — the  radius  at  the  instant  it  is  created  at  the  ocean  surface.  Lastly,  I 
assume  that  the  initial  droplet  temperature  is  Tw. 

According  to  Pruppacher  and  Klett  (1978,  p.  87) 


Ps  —  Pvv 
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It’s 

mw 
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where  mw  is  the  mass  of  pure  water  in  the  droplet  and  va  is  the  apparent  molal  volume  of 
the  solution  droplet.  I  compute  the  density  of  pure  water  in  kilograms  per  cubic  meter  at 
droplet  temperature  T  from  (Pruppacher  and  Klett  1978,  p.  86) 

p  _  999.8396  +  18.224944  7~-  7.922210x  10~3  T2  for0^r<10oC  (13a) 

1  +  1.8159725x  10  '  T 


pw  =  999.84  +  8.60  x  10'2  T  -  1.08  x  10~2  72  for  -  50  <7  <  0°C.  (13b) 
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The  apparent  molal  volume  va  in  eq  12  reflects  the  effects  of  the  dissolved  salt  on  the 
density  of  the  aqueous  solution  droplet  and  thereby  parameterizes  deviations  from  an 
ideal  solution.  Millero  (1972)  discussed  the  Masson  relationship 

=  v»o  +  V  c'/2  (14) 

which  gives  va  as  a  function  of  the  apparent  molal  volume  at  infinite  dilution  va(y  an  ex¬ 
perimental  slope  Sv*,  and  the  molar  concentration  of  the  salt  in  moles  per  liter 

c  =  10~3  .  (15) 

4  jt  r3/ 3 

Here  the  denominator  is  the  droplet  volume;  the  factor  10'3  converts  moles  per  cubic  me¬ 
ter  to  moles  per  liter  required  in  eq  14.  va0  and  Sv*  depend  on  the  solute  and  on  the  drop¬ 
let  temperature.  To  determine  what  they  are  for  NaCl,  I  plotted  Millero's  data  and  fitted 
them  with  polynomials  (Fig.  2  and  3).  In  tabulating  his  va0  data,  Millero  indicated  the  val¬ 
ues  that  he  felt  were  most  reliable.  I  fitted  the  polynomial  to  these  data  only;  the  resulting 
prediction  for  va0  in  cubic  meters  per  mole  for  0  <  T  <  55°C  is 

106ya0  =  12.97  +  0.2340  T  -  4.210  x  10~3  T2  +  2.857  x  10'5  T3.  (16) 

For  Sv*  in  m3  L1/2  mor3/2,  my  polynomial  fit  for  0  <  T  <  100°C  is 

106SV*  =  2.982  -  4.970  x  10'2  T  +  6.032  x  10"4  T2.  (17) 


Figure  2.  Apparent  molal  volume  at  infi¬ 
nite  dilution  for  NaCl  dissolved  in  pure 
water.  The  data  are  from  Millero  (1 972); 
the  curve  is  eq  16. 
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Figure  3.  Slope  in  the  Masson  equation  (eq  14)  for  NaCI  dissolved  in 
pure  ivater.  The  data  are  from  Millero  (1 972);  the  curve  is  eq  1 7. 


Notice  that  in  Figures  2  and  3  the  available  data  are  all  above  0°C;  my  polynomials,  eq 
16  and  17,  thus  have  the  same  lower  limit.  At  high  latitudes,  however,  supercooled  sea 
spray  droplets  commonly  exist.  My  only  alternative  for  modeling  these  droplets  is  to 
hope  that  eq  16  and  17  are  adequate  on  extrapolation  below  0°C. 

The  instantaneous  salinity  of  the  droplet  Sd  is  given  by 

Sd  =  -3—  .  (18) 

'"»•  +  •"* 


Consequently, 

S  = — 3 —  (19) 

'»wO+'»s 

where  ;hw0  is  the  initial  mass  of  pure  water  in  the  droplet.  Because 

(20) 

"'wO  1-S 

we  can  compute  the  initial  droplet  density  from  eq  12  if  we  approximate 

|  Kpw^-3-  (2D 

3  "'wo 

for  use  in  eq  15.  Calculating  the  initial  droplet  density  and  chosing  an  initial  radius  r0, 
we  can  find  the  initial  droplet  mass 

'"do  =  |  *  Pso  ri  ■  (22) 


'"do  =  '"wO  +  ”'s 
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we  see  from  eq  20  that  the  mass  of  salt  in  the  droplet  is 


»'s  =  s  "'do 


(24) 


which  is  assumed  to  be  constant  for  the  lifetime  ot  the  droplet. 

The  practical  osmotic  coefficient  in  eq  7  quantifies  how  dissolved  salts  affect  the 
saturation  vapor  pressure  s  at  a  plane  solution  surface.  <t> ^  is  related  to  the  water  activi¬ 
ty  aK  =  -■/‘■’•at  bv  (Pruppacher  and  Klett  1978,  p.  82) 


-  lrv?u  =  v  in  Mu.  Os . 


(25) 


Here  in  is  a  concentration  unit,  the  molality — the  number  of  moles  of  salt  dissolved  in  1 
kg  of  water 


in 


Hls 

M s  »h. 


(26) 


Low  (1969a,  1969b)  tabulated  aw  as  a  function  of  m  for  sodium  chloride  dissolved  in 
water  at  25°C.  I  converted  his  aK  values  to  d>.  according  to  eq  25  and  plotted  these  in  Fig¬ 
ure  4.  The  least-squares  polynomial  I  have  fitted  to  these  data  is 


<bs  =  0.9270  -  2.1 64  x  1  O'2  in  +  3.486  x  10'2  nr 

-  5.95b x  10"3  w/3  +  3.91  lx  10-4  w4 


(27) 


for  0  <  in  <  6. 


Figure  4.  Practical  osmotic  coefficient  of  NaCI 
dissolved  in  water  as  a  function  of  the  molali¬ 
ty.  Data  are  from  Low  (1969a,  1969b);  the 
curve  is  eq  27. 
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In  the  transfer  of  heat  and  moisture  in  the  immediate  vicinity  of  microscopic  spray 
droplets,  the  air  need  not  act  as  a  continuum;  eq  10,  therefore,  contains  modified  values 
of  the  molecular  diffusivity  of  water  vapor  in  air  DJ  and  the  thermal  conductivity  of  air 
ka.  From  Pruppacher  and  Klett  (1978,  p.  415) 


(28) 


Here  ac  =  0.036  and  A^v  =  1.3A,  =  8  x  10"8  m  are  empirical  constants,  where  X  is  the  mean 
free  path  of  air  molecules.  I  compute  the  bulk  diffusivity  Dw  at  the  droplet  temperature  T 
and  the  atmospheric  pressure  P  from  (Pruppacher  and  Klett  1978,  p.  413) 


(29) 


which  gives  Dw  in  square  meters  per  second  for  T0  =  273.15  K  and  P0  =  1013.25  hPa. 

Pruppacher  and  Klett  (1978,  p.  418)  gave  an  analogous  expression  for  the  modified 
thermal  conductivity 


_ ka  CP _ 

r  ,  (2*Ma\1/2 

r+AT  '•aTpacpa\  R  r  ' 


(30) 


Again  =  0.7  and  AT  =  2.16  x  10'7  m  are  empirical  constants.  Also  cpa  = 

1.006  x  103  J  kg'1  °C-1  is  the  specific  heat  of  air  at  constant  pressure,  and 

p, .  UH3(&)(£)  on 

is  the  density  of  air  in  kilograms  per  cubic  meter  at  the  droplet  temperature. 

I  compute  the  thermal  conductivity  of  bulk  air  ka  at  the  droplet  temperature  T  using  a 
polynomial  I  fitted  to  data  given  by  Hilsenrath  et  al.  (1960) 

ka  =  2.411  x  10'2  (1  +  3.309  x  10‘3  T-  1.441  x  10”6  T2).  (32) 

This  gives  ka  in  watts  per  meter  per  degree  Celsius  for  air  temperatures  between  -193  and 
277°C. 

Figure  5  shows  that  the  noncontinuum  effects  are,  indeed,  important,  especially  to  the 
diffusion  of  water  vapor.  In  the  figure  I  have  plotted  Dw'  and  ka  as  functions  of  the  drop¬ 
let  radius.  Even  for  a  radius  of  100  fim,  DJ  is  still  5%  below  the  bulk  diffusivity.  The  ef¬ 
fect  is  less  pronounced  on  the  diffusion  of  heat;  ka  is  the  same  as  the  bulk  thermal  con¬ 
ductivity  for  droplet  radii  of  5  fx m  and  greater. 
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Figure  5.  Modified  water  vapor  diffusivity  in  air  (eq  28)  and 
modified  thermal  conductivity  of  air  (eq  30)  as  functions  of 
the  droplet  radius.  Conditions  areT  =  0°C  and  P  =  1000  hPa. 


To  find  the  surface  tension  of  saline  water  cs  in  eq  7, 1  use  (Pruppacher  and  Klett  1978, 
p.  107) 


os  =  ow  +  2.77  x  10'5  (ms  /  mw)  (33) 

where  aw  is  the  surface  tension  of  pure  water.  According  to  Pruppacher  and  Klett  (1978, 
p.  104) 


ow  =  7.610  x  10“2  - 1.55  x  1C4  T  (34) 

which  gives  the  surface  tension  in  joules  per  square  meter  for  a  droplet  temperature  T  be¬ 
tween  -40  and  40°C. 

For  the  latent  heat  of  vaporization  of  water  Lv,  which  appears  in  eq  10  and  in  equa¬ 
tions  in  the  next  section,  I  use  (Fleagle  and  Businger  1980,  p.  113) 

Lv  =  (25.00  -  0.02274  T)xl05.  (35) 

This  gives  Lv  in  joules  per  kilogram  for  the  droplet  temperature  T  in  degrees  Celsius. 

Finally,  to  compute  the  saturation  vapor  pressure  eMt(TJ  at  air  temperature  Ta,  I  use 
Buck's  (1981)  result.  For  saturation  over  water  at  temperature  T  (in  degrees  Celsius),  he 
gave 


(T)  =  (1.0007  +  3.46x  10-6  P) 6.1121  exp  JZ502T  \  .  (36) 

r  '240.97  +  T> 

Here  the  vapor  pressure  and  the  atmospheric  pressure  P  are  both  in  hectopascals  (which 
is  equivalent  to  the  non-SI  unit  millibars). 
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THERMAL  EVOLUTION  EQUATION 


Also  in  their  Chapter  13,  Pruppacher  and  Klett  (1978)  developed  a  rate  equation  for 
the  change  in  temperature  of  an  aqueous  solution  droplet.  Their  result,  eq  (13-64),  is 

f  (Ta-T)  =  -^L_Oa'(Ta-D  +  LVDW'( Pv-Pvr)]  •  (37) 

dt  r  Pstps 

Here  c  is  the  specific  heat  of  the  spray  droplet  at  constant  pressure,  and  pv  and  pvr  are 
the  ambient  vapor  density  and  the  vapor  density  at  the  surface  of  the  droplet  respective¬ 
ly- 

Pruppacher  and  Klett  (1978)  assumed  that  the  temperature  within  a  droplet  was  uni¬ 
formly  T;  this  is  an  accurate  assumption.  We  can  show  this  by  assuming  that  there  is  no 
fluid  motion  within  a  droplet  and  then  using  the  solution  given  by  Carslaw  and  Jaeger 
(1971,  p.  233-234)  for  heat  flow  in  a  solid  sphere.  They  found  that  the  time  required  for  a 
sphere  of  radius  r  to  be  everywhere  within  4%  of  an  impulsively  applied  surface  temper¬ 
ature  is 


Dh 

where,  in  our  case,  Dh  =  1.4  x  10-7  m2/s  is  the  diffusivity  of  heat  in  sea  water  at  0°C 
(Home  1969,  p.  56).  For  1-pm  droplets,  th  =  2.9  x  10-6  s;  for  100-pm  droplets,  Th  =  2.9  x 
10“2  s.  Later,  when  I  compute  the  time  constants  xT  that  characterize  the  time  required  for 
droplets  to  come  to  thermal  equilibrium  in  air  with  temperature  Ta,  we  will  see  that  tT  is 
always  about  20  times  greater  than  Th.  Thus,  spray  droplets  are  always  well  mixed  with 
respect  to  temperature. 

The  ambient  vapor  density  in  eq  37  derives  easily  from  the  ideal  gas  law 


pv  =  15^l£s oHsLf  (39) 

R  Ta 

where  /is  again  the  relative  humidity.  Here  the  100  is  necessary  to  yield  pv  in  kilograms 
per  cubic  meter  when  is  in  hectopascals.  I  assume  that  the  air  at  the  droplet  surface  is 
in  vapor  equilibrium  with  the  droplet;  therefore  (Pruppacher  and  Klett  1978,  p.  141) 


Pvr 


lOOMwCgat  (71 
RT 


exp(y) 


where  exp(y)  gives  the  relative  humidity  at  the  droplet  surface. 


(40) 


NUMERICAL  METHODS 


The  purpose  of  the  last  two  sections  was  to  develop  equations  with  which  I  could 
model  the  size  and  thermal  evolution  of  sea  spray  droplets.  The  resulting  equations,  eq 
10  and  37,  both  have  the  form 


(41) 
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where  x  is  either  r  or  T,  and  fx  is  a  function  of  r,  T  and  environmental  parameters  such  as 
Ta,  /,  S  and  P. 

With  eq  41,  a  possible  way  to  model  the  evolution  of  x  is  with  a  single-step  finite  dif¬ 
ference  scheme  such  as  Euler's  method.  The  recursion  relation  is  (Wendroff  1969) 


*n+l  =  *n  +  Af  <t>, 


(42) 


where  Af  =  time  step 

.vn  =  current  value  of  x  (i.e.,  at  time  fn) 

•rn+i  =  predicted  value  at  the  next  time  step 

*t*n  —  /xn‘ 

On  trying  this  method  with  eq  10  and  37,  however,  I  found  it  to  diverge  before  yielding 
equilibrium  values  of  r  and  T. 

I  therefore  went  to  a  multi-step  Euler  method  with  a  predictor-corrector  technique 
(Wendroff  1969)  that  uses  both  <j>n  and  <}>n+1  to  compute  xn+r  The  series  of  recursion  rela¬ 
tions  starts  with  eq  42.  After  finding  this  initial  estimate  of  xn+1,  we  compute  (|>n+,  and 
then 


<1>  =  J^n+^n+l)  • 

We  next  recompute  xn+1  from 


(43) 


*n+i  -Xn  +  At  <fi  (44) 

and  compare  this  value  with  the  previously  computed  value  to  check  for  convergence. 
My  convergence  test  is  to  see  whether  I  [xn+, (current)  -  xn+] (previous)] /xn+1  (previous)  !< 
0.0003.  If  it  is  not,  I  use  xn+1  (current)  to  recompute  <|>n+1  and  repeat  the  process  starting 
with  eq  43. 

A  complication  of  my  analysis  is  that  I  must  model  the  evolution  of  both  r  and  T  si¬ 
multaneously.  Thus,  I  really  have  two  sets  of  identical  finite-difference  equations  like  eq 
41-44  that  I  must  handle  together— one  for  T  and  one  for  r.  The  only  essential  problem 
this  causes,  however,  is  that  the  droplet  temperature  always  reaches  its  equilibrium  value 
(T  )  long  before  the  droplet  radius  does.  At  thermal  equilibrium,  eq  37  shows  that 

T«,  =Ta  +k^wl(pv-pvr)  .  (45) 

*a 

Clearly,  since  pvr  is  still  changing,  T  changes  slowly,  though  the  droplet  is  in  virtual  ther¬ 
mal  equilibrium  with  its  environment.  To  model  these  equilibrium  temperature  changes, 
I  compute  the  droplet  temperature  at  time  step  n+ 1  from 


[  n+1 


=  1 
2 


L  D  ' 

rr>  ,  ,  ^v,n  ^w,n  /  _  >> 

I  :  'Pv  Pvr.rv 


where  the  averaging  precludes  numerical  instabilities. 


(46) 
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The  rate  at  which  spray  droplets  evolve  depends  on  their  initial  radius.  For  efficient 
computing  I  thus  chose  the  time  step  At  according  to  the  initial  radius.  By  trial  and  error  I 
found  that  a  good  initial  time  step  is 


At!  =  10ln^log<2xlo5r°^  (47) 

where  r0  is  the  initial  radius  in  meters,  and  Int{]  means  the  integer  part  of  the  expression 
in  brackets.  This  choice  does  have  some  physical  basis,  because  the  rate  at  which  a  drop¬ 
let  exhanges  heat  and  moisture  is  almost  directly  proportional  to  its  surface  area,  4rcr.  I 
computed  the  running  time  from 

'n+l  =  'n  +  Afn.  (48) 

Time  steps  after  the  first  one  are  logarithmic;  I  set  them  according  to 

A  fn  =  0.1  fn .  (49) 

To  test  the  equations  that  Pruppacher  and  Klett  (1978)  developed,  my  evaluations  of 
the  physical  and  chemical  variables  involved,  and  the  numerical  method,  I  compared  the 
data  of  El  Golli  et  al.  (1974)  and  Ranz  and  Marshall  (1952)  with  my  model  predictions 
(Fig.  6).  El  Golli  et  al.  observed  the  evaporation  of  saline  droplets  as  they  were  advected 
by  an  air  stream  through  a  pipe  with  controlled  temperature  and  humidity.  Although  the 
droplets  that  El  Golli  et  al.  observed  had  initial  radii  near  the  median  of  sea  spray  drop¬ 
lets — about  8  /an — their  experimental  parameters  are  not  especially  representative  of  a 
sea  spray  environment.  Typically  for  sea  spray,  S  -  34  %o  and  RH  >  80%.  El  Golli  et  al. 
observed  droplets  in  an  environment  with  RH  =  80%,  but  for  these  the  initial  salinity  was 
only  2  %o  (Fig.  6c).  In  their  highest  salinity  case,  S  =  29  %o,  the  relative  humidity  was  only 
29%  (Fig.  6a).  Fortunately,  nothing  in  the  equations  that  I  presented  in  the  Size  Evolution 
Equation  and  Thermal  Evolution  Equation  sections  makes  them  exclusively  applicable  to  the 
high  salinity  sea  spray  environment;  they  should  also  be  accurate  at  lower  salinity. 

We  must,  however,  be  careful  treating  low-humidity  cases.  At  low  humidity  a  droplet 
may  lose  enough  water  by  evaporation  to  become  saturated  with  salt.  At  0°C,  1  L  of  pure 
water  can  dissolve  0.357  kg  of  sodium  chloride;  thus,  the  maximum  molality  of  a  saline 
solution  is  6.11  (see  eq  26).  Such  a  saturated  solution  droplet  is  in  moisture  equilibrium 
with  air  having  a  relative  humidity  of  about  75%  (Twomey  1954).  At  lower  humidity,  the 
salt  will  eventually  crystallize  and  the  water  will  all  evaporate.  For  the  two  low-humidity 
data  sets  collected  by  El  Golli  et  al.  (1974)  (Fig.  6a  and  b),  I  therefore  stopped  the  calcula¬ 
tions  when  enough  water  had  evaporated  to  increase  the  droplet  molality  to  6.11. 

One  ambiguity  that  I  found  in  the  experimental  data  of  El  Golli  et  al.  (1974)  is  that 
they  reported  only  one  temperature  for  each  set  of  experimental  conditions.  They  said 
this  was  an  "average"  temperature  but  did  not  clarify  whether  it  was  the  average  air  tem¬ 
perature  in  their  device,  the  average  droplet  temperature  during  the  experimental  run,  or 
some  other  average.  In  their  experiments  the  droplets  were  in  thermal  equilibrium  by  the 
time  they  started  observing  their  evaporation.  From  eq  45  we  see  that,  because  of  the 
evaporation,  the  droplets  could  not  have  been  at  Ta — they  had  to  be  cooler.  A  second 
problem  with  their  reporting  of  the  experimental  temperature  is  that  the  moisture  in  the 
air  at  a  given  relative  humidity  is  a  strong  function  of  the  air  temperature.  Consequently, 
since  El  Golli  et  al.  reported  RH  and  some  "average"  temperature,  the  comparison  de¬ 
pends  critically  on  what  I  choose  for  Ta.  The  conditions  given  in  Figure  6  show  the  tem- 
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Time  (s) 


a.  Data  from  Table  2  in  El  Golli  et  al.  (1974).  Model  con¬ 
ditions  are  RH  =  29%,  S  =  29  %c,  Tcq  =  26°C,  T „  =  39  °C 
and  P  =  1000  hPa. 


b.  Data  from  Table  3  in  El  Golli  et  al.  (1974).  Model  con¬ 
ditions  are  RH  =  32%,  S  =  10  %c,  Tcq  =  28°C,  T„  =  41  °C, 
and  P  =  1000  hPa. 

Figure  6.  Model  computations  compared  -with  experiment¬ 
al  data. 
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peratures  that  yield  the  best  agreement  with  the  data  from  El  Golli  et  al.  It  turned  out  that 
in  my  notation  their  reported  "average"  temperature  is  T  ,  the  droplet  temperature  at 
equilibrium.  I,  thus,  chose  Ta  to  yield  this  Tecj  value,  given  the  RH  and  S  values  reported 
by  El  Golli  et  al. 

Aside  from  Figure  6c,  where  my  model  overestimates  the  experimental  droplet  size  at 
any  given  time  by  about  30%,  on  average,  my  model  computations  agree  quite  well  with 
the  size  evolution  data  reported  by  El  Golli  et  al.  (1974). 

Data  presented  graphically  by  Ranz  and  Marshall  (1952,  their  Fig.  15b)  provide  anoth¬ 
er  important  test  of  my  model.  Although  they  considered  large  and  unnaturally  saline 
droplets  (S  =  214  %o),  my  model  simulates  their  data  well.  Figure  6d  shows  the  compari¬ 
son;  my  model  is  always  within  9%  of  their  experimental  data.  Ranz  and  Marshall  did 
not  specify  the  ambient  relative  humidity  for  their  experiment  but  assumed  it  to  be  near 
zero.  By  trial  and  error,  I  found  that  a  relative  humidity  of  8%  produced  the  best  agree¬ 
ment  between  the  model  values  of  7^  and  the  droplet  temperatures  measured  by  Ranz 
and  Marshall. 

My  analysis  of  the  Ranz  and  Marshall  (1952)  data  also  substantiates  some  of  the  un¬ 
derlying  physics  in  my  model.  Cheng  et  al.  (1988)  recently  reported  producing  hollow, 
spherical  salt  particles  by  evaporating  sea  water  droplets  in  a  laboratory  and  noted  that 
Ranz  and  Marshall  observed  similar  evaporative  particles.  Cheng  et  al.  speculated  that 
the  hollow  spheres  might  form  because  sea  salts  have  low  diffusivity  in  water  and,  thus, 
might  concentrate  at  the  surface  of  evaporating  droplets.  My  model  would  be  in  trouble 
if  this  mechanism  is  at  work  because  it  assumes  that  the  salt  within  the  droplet  is  well 
mixed.  Model  simulations  of  the  Ranz-Marshall  data,  fortunately,  repudiate  this  mecha¬ 
nism  and,  thus,  confirm  my  assumption  that  the  salt  is  well  mixed.  In  their  Figure  15b, 
Ranz  and  Marshall  noted  that  the  first  salt  crystals  began  forming  near  a  time  of  380  s.  In 
my  model  simulations  of  their  data,  the  droplet  molality  reached  6.11 — the  salt  concen¬ 
tration  at  saturation — near  390  s.  (In  Figure  6d  I  plot  no  data  or  model  results  beyond  this 
limit.)  This  excellent  agreement  supports  my  assumption  that  the  salt  in  the  droplets  is 
well  mixed.  If  salt  concentrates  near  the  surface  of  evaporating  droplets,  as  Cheng  et  al. 
hypothesized,  Ranz  and  Marshall  should  have  observed  crystals  forming  much  earlier. 

Using  an  analysis  similar  to  that  given  in  the  Thermal  Evolution  Equation  section,  I  can 
look  closer  at  the  hypothesis  by  Cheng  et  al.  (1988)  that  salt  diffusion  in  evaporating 
droplets  is  so  slow  that  salt  concentrates  near  the  surface.  In  the  Thermal  Evolution  Equa¬ 
tion  section,  I  presented  Carslaw  and  Jaeger's  (1971,  p.  233-234)  solution  for  heat  diffu¬ 
sion  in  a  solid  sphere.  Salt  diffusion  is  an  analogous  problem.  An  incremental  surface 
evaporation  is  equivalent  to  an  incremental  injection  of  salt  at  the  droplet  surface.  This 
salt  will  tend  to  diffuse  throughout  the  droplet  by  molecular  processes.  Unlike  the  heat 
diffusion  problem  in  which  the  surface  was  clamped  at  the  new  temperature,  however, 
the  salt  will  become  diluted  as  it  diffuses  throughout  the  droplet;  the  final  concentration 
will,  therefore,  be  everywhere  less  than  the  original  surface  concentration.  Consequently, 
a  time  constant  like  eq  38  will  be  an  upper  bound  on  the  time  required  (ts)  for  a  surface 
salt  injection  to  diffuse  throughout  the  droplet 


<  0-4  d 
Ds 


(50) 


Here  D s  =  6.8  x  10  10  m2/s  is  the  molecular  diffusivity  of  NaCl  in  sea  water  at  0°C  (Horne 
1969,  p.  56).  For  1-jim  droplets,  ts  <  5.9  x  10"*  s;  for  100-pm  droplets,  <  5.9  s.  We  will  see 
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from  computations  that  I  will  discuss  later  that  these  diffusion  times  are  always  two  or¬ 
ders  of  magnitude  less  that  the  time  constant  tf,  which  characterizes  the  time  required  for 
sea  spray  droplets  to  reach  moisture  equilibrium  with  their  environment.  In  other  words, 
salt  diffusion  within  the  droplets  always  seems  to  be  fast  enough  to  maintain  a  well- 
mixed  interior  during  sea  spray  evaporation  in  a  typical  marine  environment. 

With  the  success  of  my  model  in  simulating  the  data  from  El  Golli  et  al.  (1974)  and  the 
data  and  observations  of  Ranz  and  Marshall  (1952),  1  feel  that  we  can  use  it  with  confi¬ 
dence  to  investigate  the  thermal  and  size  evolution  of  sea  spray  droplets. 


THE  GENERATION  OF  SEA  SPRAY 

As  I  hinted  in  the  Introduction,  sea  spray  droplets  form  by  three  distinct  mechanisms. 
When  an  air  bubble  rises  to  the  ocean  surface,  the  bubole  cap  thins  and  finally  ruptures, 
throwing  micrometer  and  submicrometer  sized  film  droplets  into  the  air.  As  the  bubble 
cavity  then  collapses,  one  to  several  jet  droplets  shoot  into  the  air  from  its  base.  The  top  jet 
droplet  is  about  one-tenth  the  radius  of  the  original  bubble;  lower  jet  droplets  are  gener¬ 
ally  larger  (Blanchard  and  Syzdek  1988).  Blanchard  (1963, 1983),  Cipriano  and  Blanchard 
(1981)  and  Woolf  et  al.  (1987)  discussed  these  mechanisms;  Kientzler  et  al.  (1954),  Day 
(1964)  and  Resch  et  al.  (1986)  showed  dramatic  photos  of  bursting  bubbles  and  the  result¬ 
ing  spray.  When  the  10-m  wind  reaches  speeds  higher  than  about  9  m/s,  it  has  the  ener¬ 
gy  to  tear  off  the  wave  crest  and  thereby  produce  spume  droplets  (Monahan  et  al.  1983). 
These  are  typically  the  largest  sea  spray  droplets,  with  radii  generally  greater  that  10  pm. 

Information  on  the  generation  rate  of  sea  spray  droplets  in  a  given  size  range  comes 
from  many  sources.  Monahan  (1968),  Preobrazhenskii  (1973),  Monahan  et  al.  (1983)  and 
Wu  et  al.  (1984)  reported  measurements  of  the  spray  droplet  size  distribution  at  the  sea 
surface.  Wu  (1973),  Lai  and  Shemdin  (1974),  Wang  and  Street  (1978)  and  Mestayer  and 
Lefauconnier  (1988)  measured  the  spray  size  distribution  in  wind-wave  tunnels.  Mona¬ 
han  et  al.  (1982)  measured  it  in  a  whitecap  simulation  tank. 

Monahan  et  al.  (1986)  synthesized  much  of  these  data  in  a  single  spray  generation 
function  by  recognizing  that  two  mechanical  processes  produce  sea  spray — bursting  bub¬ 
bles  (film  and  jet  droplets)  and  tearing  of  wave  crests  by  the  wind  (spume  droplets). 
Their  spray  generation  function  dF/drm  gives  the  number  of  spray  droplets  produced  per 
unit  surface  area  per  second  per  micrometer  increment  in  droplet  radius;  the  units  are 
thus  m  2  s'1  /im"1  (where  number  is  implicit  in  the  numerator).  In  dF/drm,  r80  denotes  that 
all  the  observational  data  were  converted  to  an  ambient  relative  humidity  of  80%.  Mona¬ 
han  et  al.  obtained  dF/drm  from  the  sum  of  two  terms 

dF_  =  dFo+dFL  (51) 

^r80  ^r80  ^r80 


where  dF0/drm  and  dFl/drm  are  the  contributions  from  bubble  and  spume  droplets  re¬ 
spectively. 

For  dFJdr^  Monahan  et  al.  (1986)  gave 

=  1.3731/ 341  r'^(l  +  0.057r^°5)l0U9l'xp('B2>  (52a) 
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with 


B  ~0~650  (0'380  ~  logrso)  '  (52b) 

Here  U  is  the  10-m  wind  speed  in  meters  per  second. 

For  dFx/drS0  Monahan  et  al.  (1986)  gave 


4L >  =o 

dr  so 

for  r80  <  10^m 

(53a) 

=  8.60  x  10-6  exp  (2.08  U)  r ^ 

for  10  S  rjjo  <  75  jjm 

(53b) 

=  4.83  x  10"2  exp(2.08  U)  r^ 

for  75  S  rgo  <  100  /Jm 

(53c) 

=  4.83  x  10b  exp(2.08  U)  r^ 

for  100  ^im  <  rs0 

(53d) 

Both  eq  52  and  53  are  accurate  only  for  a  wind  speed  U  that  is,  nominally,  20  m/s  or  less. 

Because  my  ultimate  interests  require  that  I  track  the  evolution  of  spray  droplets  from 
the  instant  they  form,  a  spray  generation  function  in  terms  of  r80  is  inadequate.  From 
dF/drm  I  can,  however,  find  dF/dr,  the  spray  generation  function  in  terms  of  the  droplet 
radius  at  formation  r(y  from 

d£_  _  d]jo_dF_ 

dr  dr0  dr w  (54) 

Of  course,  all  the  data  that  Monahan  et  al.  (1986)  considered  were  not  obtained  at  a 
relative  humidity  of  80%;  to  have  compatible  data  sets  they  had  to  use  a  transformation 
like  eq  54  to  convert  the  various  observations  to  a  common  humidity  (they  chose  80%). 
To  find  drm/drg,  Monahan  et  al.  (see  also  Fairall  et  al.  1983)  used  Fitzgerald's  (1975)  result 

req  =  «(/)  rdp(f)  .  (55) 


Here  riH)  is  the  equilibrium  droplet  radius  at  relative  humidity  /,  rd  is  the  radius  of  the  dry 
salt  nucleus,  and  for  NaCl 


where 


gif)  =  1.058 


for/  <097 


(57a) 


=  1.058 


0.0155  if  -  0.97) 
1.02  -f'A 


for  0.97  <  /  <  0.995 


(57b) 
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and 

(3(f)  =  exp(«)  for  0.80  </<  0.995 .  (58) 

For  eq  58,  Fitzgerald  gave  the  lower  limit  of  applicability  as  /  =  0.81;  but  I  see  nothing  in 
his  work  to  preclude  us  from  extending  eq  58  down  to  /=  0.80. 

I  assume  that  when  spray  droplets  form,  they  are  injected  initially  into  an  environ¬ 
ment  that  is  saturated  with  water  vapor — the  lowest  few  centimeters  of  air  near  the  sea 
surface.  But  because  the  dissolved  salts  depress  the  saturation  vapor  pressure  compared 
to  that  over  a  plane  surface  of  pure  water,  the  initial  relative  humidity  is  not  100%.  Rath¬ 
er,  it  is  (Roll  1965,  p.  262;  Home  1969,  p.  332) 

/0  =  1  -0.000537  S  (59) 

where  S  is  the  sea  surface  salinity  in  parts  per  thousand. 

Knowing  this  initial  humidity,  we  can  use  eq  55  to  relate  r80  to  r0  and,  thus,  find 
drm/drQ.  The  procedure  is  to  write  the  two  equations 

'o  =  a  (/o)  rj <ffl)  (60) 


rso  =  a  (0.8)  rd  (0'8)  (61) 

and  to  eliminate  rd.  Figure  7  shows  versus  rQ  computed  from  both  Fitzgerald's  (1975) 
equations  and  from  my  model.  A  log-log,  least-squares  fit  of  the  r80-r0  data  computed 
with  Fitzgerald's  equations  is 

r80  =  0.51 75r  o'9756  (62) 

with  both  r80  and  r0  in  micrometers.  The  least-squares  fit  to  the  results  from  my  model  is 

r80  =  0.4954r  J,0029  .  (63) 


Looking  at  these  two  equations  and  at  Figure  7,  we  see  little  difference  between  the 
two  methods  of  computing  rm.  But,  as  Figure  8  shows,  there  are  significant  differences 
between  Fitzgerald's  (1975)  relations  and  my  model.  In  Figure  8 1  plot  the  ratio  r0/r^  ver¬ 
sus  r0  for  several  values  of  the  ambient  relative  humidity.  Here  is  the  equilibrium  radi¬ 
us  at  the  specified  humidity;  for  exampie,  for  RH  =  80%,  =  r80  and  that  set  of  lines  in 

Figure  8  contains  the  same  information  as  Figure  7.  Clearly,  in  Figure  8,  r0/r ■  depends 
much  more  on  r0  in  Fitzgerald's  model  than  in  my  model.  It  is  not,  however,  clear  which 
model  is  correct.  I  can  only  say  that  I  feel  mine  is  because  it  computes  from  the  funda¬ 
mental  microphysical  equations.  Fitzgerald's  relations,  on  the  other  hand,  are  approxi¬ 
mate  fits  to  solutions  of  the  microphysical  equations  in  slightly  different  form,  though  he 
claimed  that  the  approximations  are  good  to  within  5%.  I  also  may  be  extrapolating  his 
results  too  far;  this  stated  accuracy  applies  only  when  r0  is  10-20  pm  or  less. 

Since  my  calculations  of  r0/ r^  yield  values  that  are  basically  constant  over  a  wide  ra¬ 
dius  range,  they  have  an  advantage  over  computations  based  on  Fitzgerald's  (1975)  mod- 
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Figure  7.  Comparison  between  r80  and  r0  based  on 
Fitzgerald's  (1975)  relations  and  on  the  present  model. 


Figure  8.  Ratio  computed  with  Fitzgerald's  (1975)  rela¬ 
tions  and  with  the  present  model  for  several  values  of  the  am¬ 
bient  relative  humidity. 


el.  In  fact.  Figure  8  provides  a  rationale  for  reducing  experimental  sea  spray  data  to  a 
standard  humidity  of  80%:  at  RH  =  80%  a  spray  droplet  has  a  radius  that  is  one-half  the 
droplet  radius  at  formation. 

Although  computations  of  rm  versus  r0  based  on  my  model  seem  aesthetically  better 
than  those  based  on  Fitzgerald's  (1975)  model,  because  Monahan  et  al.  (1986)  used  Fitz¬ 
gerald's  relations  to  reduce  their  experimental  data  to  a  common  humidity,  I  defer  to  the 
Fitzgerald  model  in  computing  drm/ dr0  for  use  in  eq  54.  Therefore,  from  eq  62 

drso  -  fl  S049r'0  0244  ,,,, 

— - u.ou4*r0  (64) 

ur  o 

Using  this  and  eq  51-54, 1  can  compute  the  rate  at  which  spray  droplets  are  formed  as 
a  function  of  the  wind  speed.  Notice,  none  of  the  equations  I  use  for  this  calculation  de¬ 
pend  on  the  humidity,  salinity,  or  air  or  water  temperatures.  Some  preliminary  evidence 
does,  however,  suggest  that  the  spray  generation  function  depends  on  the  water  temper¬ 
ature.*  Figure  9  shows  dF/dr,  the  spray  generation  function,  as  a  function  of  r0  and  the 
wind  speed.  Not  unexpectedly,  the  number  of  droplets  produced  with  a  given  radius  in¬ 
creases  monotonically  with  the  wind  speed.  This  is  the  bursting  bubble  mechanism. 

Figure  9  also  shows  clearly  that  at  a  wind  speed  of  about  9  m/s  a  new  spray  genera¬ 
tion  mechanism  becomes  effective.  This  is  spume  generation.  The  wind  is  now  powerful 

enough  to  tear  off  the  wave  crests  and  to  produce  many  large  spray  droplets. 

The  basic  purpose  of  this  section  was  to  identify  the  size  range  of  spray  droplets  with 
which  we  must  be  concerned.  Figure  9  delineates  that  size  range.  Although  the  spray 
generation  function  increases  monotonically  with  decreasing  radius  in  Figure  9,  we  can 
call  r0  =  0.5  jum  a  lower  limit,  because  the  mass  carried  by  these  smallest  droplets  be¬ 
comes  negligible.  Remember,  the  mass  of  a  spray  droplet  goes  as  r03.  A  reasonable  upper 
size  limit  is  r0  =  500  fim.  As  I  will  show  later,  500-/im  spray  droplets  remain  in  the  air  less 
than  a  second  before  falling  back  into  the  sea.  Henceforth,  I  will  limit  my  discussion  of 
sea  spray  to  droplets  with  an  initial  radius  between  0.5  and  500  jum. 


EXAMPLES  OF  THERMAL  AND  SIZE  EVOLUTION 

As  I  explained  in  the  Introduction,  I  am  interested  in  spray  generation  associated  with 
polar  lows.  Since  these  occur  at  high  northern  latitudes,  the  sea  surface  temperature  Tw  is 
typically  0°C.  If  the  wind  is  blowing  off  the  pack  ice,  the  air  temperature  may  be  as  low 
as  -20°C.  For  southerly  winds,  the  air  is  warmer;  it  could  reach,  say,  5°C.  The  surface  sa¬ 
linity  S  of  the  high  latitude  ocean  is  commonly  about  34  %o.  These  are  the  environmental 
parameters  that  I  will  use  for  my  studies  of  the  thermal  and  size  evolution  of  sea  spray 
droplets. 

Thermal  evolution 

Figure  10  shows  three  plots  of  the  evolution  of  the  droplet  temperature  T  with  time.  I 
always  assume  that  a  spray  droplet  initially  has  the  temperature  of  the  surface  water  Tw. 

is  the  equilibrium  temperature  of  the  droplet  in  the  specified  environment  and  is  giv¬ 
en  by  eq  45.  In  Figure  10  I  plot  the  nondimensional  temperature  [  7(0-7^} /[T^-T^]  as  a 


*  Personal  communication  with  E.  C.  Monahan,  University  of  Connecticut,  1987. 
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Figure  9.  Spray  generation  function  (eq  54)  for  various  'values  of 
the  10-in  wind  speed. 


function  of  the  nondimensional  time  f/xT.  Here  xT  is  a  time  constant  equal  to  the  time  re¬ 
quired  for  the  spray  droplet  to  reach  ir'  of  its  equilibrium  temperature.  That  is 

Tk  -  T0l| 

In  the  next  section  I  will  quantify  xT. 

The  thermal  evolution  depicted  in  Figure  10  is  virtually  independent  of  the  relative 
humidity.  The  spray  droplets  reach  thermal  equilibrium  so  quickly  that  no  mass  has  yet 
been  exchanged. 

The  three  droplets  modeled  in  Figure  10  have  nearly  identical  thermal  evolution 
curves.  Thus,  the  basic  character  of  thermal  transfer  from  a  sea  spray  droplet  does  not  de- 


(65) 
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i/rT  i/tt 


b.  r0=  10  \xm,Ta  =  -20°C.  c.  r„  =  10  \im,  Tn  =  S'C. 

Figure  10.  Model  calculations  of  the  thermal  evolution  of  spray  droplets.  Environmental  condi¬ 
tions  arc  T.,,  =  D°C,S  =  34  %c  and  P  =  1 000  hPa. 


pend  strongly  on  its  size,  on  its  initial  temperature  or  on  .the  ambient  temperature  and 
humidity.  It  does  not  even  depend  strongly  on  the  direction  of  the  transfer.  In  Figure  10c 
the  droplet  is  warming,  and  the  nondimensional  temperature  generally  falls  below  the 
line 


T(t)-T 


eq 


TW~T , 


eq 


exp  (— f /xT)  • 


(66) 


Although  the  nondimensional  temperature  is  above  this  line  for  the  two  cooling  cases 
(Fig.  10a  and  10b),  the  deviations  from  exp(-f/tT)  are  not  large.  Consequently,  if  we 
know  TeC[ — which  we  do  from  eq  45 — and  tT,  exp(-f/xT)  is  a  fairly  accurate  model  for  the 
thermal  evolution  of  spray  droplets,  regardless  of  the  environmental  conditions  (see  also 
Pruppacher  and  Klett  1978,  p.  447). 

Size  evolution 

Figure  11  shows  spray  droplet  size  evolution  for  the  three  cases  shown  in  Figure  10.  In 
analogy  with  the  temperature,  I  plot  the  nondimensional  radius  [r(f)-req]/[r0-re[]]  versus 
the  nondimensional  time  t/xr  Here  tr  is  a  time  constant  equal  to  the  time  required  for  the 
droplet  to  come  to  within  e~]  of  its  equilibrium  radius.  That  is 


r(xr)  -  rcq 
ro-rcq 


(67) 


I  will  discuss  values  of  xx  in  the  next  section. 

Figure  11  shows  that  the  size  evolution  depends  strongly  on  the  relative  humidity. 
The  lower  the  relative  humidity,  the  more  rapid  the  moisture  exchange  at  the  droplet  sur¬ 
face. 

Clearly,  the  model 

— — ^  =  exp  (~t/xr)  (68) 

'o-rc  q 


is  not  as  good  for  size  as  it  is  for  droplet  temperature.  Only  for  t/xx  <  1  is  eq  68  a  reason¬ 
able  model  of  the  droplet  radius. 

Figure  12  presents  the  same  results  shown  in  Figure  11  but  in  a  form  that  emphasizes 
how  much  a  spray  droplet  actually  changes  in  size  and  how  r  depends  on  the  relative 
humidity.  The  figure,  therefore,  reiterates  some  of  the  information  contained  in  Figure  8: 
for  RH  =  80%,  r0/ =  2.0;  for  RH  =  90%,  /'0/ret)  =  1 .7;  and  for  RH  =  95%,  /-(,/req  =  1 .3. 

On  comparing  Figures  11  and  12,  we  see  that  r{t)/r0  is  a  more  transparent  way  to  non- 
dimensionalize  size  evolution  than  is  [r(f)-r  ]/[r0-r  ];  the  three  plots  in  Figure  12  are 
nearly  identical  despite  the  diverse  conditions  that  they  represent.  In  Figure  11  the  three 
plots  are  noticeably  different.  Evidently,  when  we  parameterize  size  evolution  as  r(f)/r() 
versus  t/xr,  the  initial  droplet  size  and  the  ambient  air  and  water  temperatures  do  not 
strongly  affect  the  nature  of  the  moisture  transfer.  As  I  will  show  in  the  next  section,  they 
do,  however,  affect  the  rate  of  the  evolution — modeled  by  tr. 
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b.  r0  =  10  pm,  Ta  =  -20°C. 


c.  r0  =  10  pm,Ta  =  5°C. 


Figure  11,  Model  calculations  of  the  size  evolution  of  spray  droplets  plotted  nondimensionally  as 
lr(tF-r  ]/lr0- r^].  Environmental  conditions  are  the  same  as  for  Figure  10. 


Figure  12.  Model  computations  of  the  size  evolution  of 
spray  droplets  plotted  nondimensionally  as  r(t)fr0. 
Environmental  conditions  are  the  same  as  for  Figure  1 0. 


c.  r0=  10[xm,Tn  =  5°C. 
Figure  12  (cont'd). 


Figure  12  shows  that,  when  the  relative  humidity  is  98%,  the  spray  droplets  grow 
somewhat,  regardless  of  the  direction  of  the  thermal  transfer.  According  to  eq  10,  droplet 
growth  for  RH  =  98%  can  occur  only  when 

i /  <  0.98  -1  =  -0.02.  (69) 

From  eq  7  we  can  easily  demonstrate  that  y  is  typically  somewhat  smaller  than  -0.02  ex¬ 
cept  for  the  smallest  droplets,  for  which  surface  tension  is  an  important  effect.  Surface 
tension  exerts  its  influence  when  the  initial  droplet  radius  is  about  1  /urn  or  less.  Figure  12 
will  thus  not  give  an  accurate  picture  of  the  evolution  of  droplets  of  this  size;  when  the 
relative  humidity  is  98%,  they  will  reach  an  equilibrium  size  slightly  less  than  r0 


TIME  CONSTANTS 

The  time  constants  tT  and  tr  that  I  defined  in  the  last  section  parameterize  how  quick¬ 
ly  a  spray  droplet  comes  to  equilibrium  with  its  environment.  With  my  model  I  can  com¬ 
pute  ty  and  tr  for  any  initial  radius  and  for  any  set  of  ambient  conditions  with  RH  >  75%. 
Figure  13  shows  plots  of  Xj  and  Tr  for  droplets  with  initial  radii  from  0.5  to  500  /urn  and 
for  a  host  of  ambient  conditions.  Notice  the  tt  and  Tr  values.  A  1-qm  droplet  is  within  c~] 
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of  its  equilibrium  temperature  in  less  than  lCT*  s;  a  lOO-pm  droplet,  in  less  than  1  s.  A  1- 
/im  droplet  is  within  e~l  of  its  equilibrium  radius  in  about  0.5  s;  a  100-pm  droplet,  in 
about  1000  s.  Clearly,  the  time  constants  depend  on  the  droplet  radius  and  on  the  constit¬ 
uent  being  exchanged. 

The  Zj  line  in  each  plot  of  Figure  13  depends  negligibly  on  the  humidity  and  little  on 
the  air  or  water  temperatures.  The  figure  itself  explains  the  lack  of  humidity  dependence. 
Because  ^  is  about  three  orders  of  magnitude  less  than  tr,  moisture  exchange  at  the  drop¬ 
let  surface  can  have  no  influence  on  thermal  exchange;  a  droplet  reaches  thermal  equilib¬ 
rium  long  before  moisture  diffusion  begins  in  earnest  at  the  droplet  surface.  Conversely, 
because  a  droplet  is  in  thermal  equilibrium  with  the  air  long  before  moisture  transfer  be¬ 
gins,  the  initial  droplet-air  temperature  difference  can  have  no  effect  on  that  transfer. 
Hence,  the  choice  of  Tw  is  immaterial  to  the  size  evolution  calculations. 

In  each  plot  in  Figure  13,  the  Xj  line  is  almost  straight  on  these  log-log  plots  and  has  a 
slope  of  2.  That  is, 

tT-r02.  (70) 


a.  Ta  =  -20°C.  b.  T a  =  -10°C. 


Figure  13.  Time  constants  for  thermal  (rT)  and  moisture  (rj  exchange.  Environmental  condi¬ 
tions  are  Tw  =  0°C,  S  =  34%c  and  P  =  1000  hPa. 
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Time  (s) 


i  10  'iod ~  *iOCO 

'o  (  1 


c.  T a  =  0°C. 


d.  T a  =  5°C. 


Figure  13  {coni' A). 


Something  inherent  in  the  thermal  evolution  equation  must  be  responsible  for  this  simple 
result.  So  consider  eq  37  again.  From  eq  45  we  can  approximate  this  as 

£(Teq-T)  .  (71) 

dt  r»  PsCps 

The  Tjq  substitution  on  the  right-hand  side  here  is  an  approximation  because  Lv  and  DJ 
in  eq  37  depend  on  T.  Therefore,  eq  45  is  really  valid  only  when  all  the  variables  are  eval¬ 
uated  at  T  .  In  eq  71  I  also  substituted  T^  for  Ta  on  the  left-hand  side  since  both  are  as¬ 
sumed  to  be  independent  of  time.  Because  the  droplet  temperature  reaches  equilibrium 
long  before  any  significant  moisture  transfer  occurs,  r  and  ka\  which  are  also  time  depen¬ 
dent,  are  essentially  constants  while  a  droplet  comes  to  thermal  equilibrium.  ps,  the  solu¬ 
tion  density,  is  also  time  dependent  in  my  model;  but  since  that  time  dependence  is 
through  a  weak  dependence  on  T,  I  ignore  it,  too.  Therefore,  eq  71  integrates  easily  to 


In 


7X0 -T*, 


-3k/t 

r  0  Ps^ps 


(72) 
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When  t  =  Xj,  IT(tt)-Tn]/[Tw-T^]  =  «r>.  Therefore, 


-■-/j'Ttc  =  -1  (73) 

''OPsCps 


or 


TTc  - 


Ps  ^*ps  ro 

3  k/ 


(74) 


Here  the  subscript  c  denotes  that  this  is  a  time  constant  calculated  with  the  approxima¬ 
tion  eq  74  as  opposed  to  the  value  xT  based  on  the  full  model. 

Figure  14  shows  that  eq  74  may  be  an  adequate  approximation  for  some  applications. 
In  the  figure  I  plot  (xTc-xT)/xT,  which  is  the  relative  error  in  the  xTc  estimate.  xTc  always 
overestimates  xT,  but  the  relative  error  for  all  radii  of  interest  is  less  than  50%.  For  the 
smallest  droplets,  that  error  is  even  less.  Again,  these  errors  result  because  in  reducing  cq 
37  to  eq  74  I  made  several  approximations;  most  involved  ignoring  the  dependence  on 
the  instantaneous  droplet  temperature  of  almost  all  the  parameters  in  eq  37. 

In  each  plot  in  Figure  13,  xr  increases  as  the  relative  humidity  increases.  The  reason  is 
obvious  in  eq  10.  The  larger  I/  - 1  I  is,  the  faster  the  rate  at  which  a  droplet  radius  changes. 

I  have  ordered  the  individual  plots  in  Figure  13  such  that  13a  has  the  largest  values  of 
xr  at  a  given  radius  and  13d  has  the  smallest.  This  sequence  also  reflects  air  temperatures 
that  go  -20,  -10,  0,  5°C.  The  obvious  conclusion  is  that  xr  decreases  with  an  increase  in 
ambient  temperature.  What's  happening  physically  is  that  a  warmer  droplet  exchanges 
vapor  with  its  environment  more  rapidly  than  does  a  cooler  one  (Bohren  1986).  Mathe¬ 
matically,  the  effect  manifests  through  the  ^(T,)  term  in  eq  1  and  10. 


Figure  14.  Comparison  of  calculated  ( rT.,  from  eq  74)  and  model  ( rr) 
values  of  the  time  constant  for  thermal  equilibrium.  For  the  model  calcula¬ 
tions,  environmental  conditions  are  T(l,  =  0°C,  Ta  =  -10  °C,  RH  =  70%  (not 
important),  S  =  34  9 It  and  P  =  1000  hPa.  For  computing  rTi„  all  tempera¬ 
ture-dependent  parameters  were  evaluated  at  0°C. 
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The  largest  values  for  tr  depicted  in  Figure  13  are  many  hours.  It  is  very  unlikely  that 
the  environment  in  which  sea  spray  droplets  evolve  will  remain  constant  this  long.  And 
even  if  the  environment  remains  constant,  it  is  even  more  unlikely  that  the  largest  drop¬ 
lets  will  stay  suspended  in  the  air  long  enough  to  approach  equilibrium. 

We  can  make  a  rough  estimate  of  this  atmospheric  residence  time  by  computing  the 
Stokes  fall  speed  of  each  droplet.  For  an  atmosphere  at  rest,  the  force  balance  for  a  spray 
droplet  of  radius  r  is 


~  rc  r3  (ps  -Pa)  $  ~CD  pa  (*  r 2)  u1  =  1 71  r3  psfl 


(75) 


where  y  =  acceleration  of  gravity 

CD  =  drag  coefficient  of  the  droplet 

u  =  droplet  fall  speed 

a  =  droplet  acceleration. 

In  eq  75  the  left-most  term  is  the  gravitational  force  on  the  droplet — modified  by  buoyan¬ 
cy  effects.  The  middle  term  is  the  frictional  drag.  The  term  on  the  right  is  the  resulting  ac¬ 
celeration.  When  a  droplet  reaches  its  terminal  fall  speed  uf,  the  right-hand  term  in  eq  75 
is  zero.  Then 


u 2  _  4  Kps  -Paly 


3  CD  pa 


(76) 


Schlichting  (1968,  p.  106)  and  Batchelor  (1970,  p.  233)  gave  the  drag  coefficient  for  a 
sphere  as 


Cr,  -  6  v* 
C-D - 

u  r 


(77) 


(Notice,  my  definition  of  CD  in  eq  76  reflects  a  value  half  as  large  as  that  used  by  Schlich¬ 
ting  and  Batchelor.)  The  kinematic  viscosity  of  air  va  in  eq  77  is  a  function  of  the  air  tem¬ 
perature  Ta;  to  compute  va  for  typical  sea  level  pressures  near  1000  hPa,  I  fitted  the  data 
of  Hilsenrath  et  al.  (1960)  with  the  polynomial 


va(7)  =  1.326  xl0-3(l  +  6.542  xl0-3T 


+  8.301  x  10^  T2  -  4.840  x  10'v  D  . 


,-y  -r-b 


(78) 


This  yields  va  in  square  meters  per  second  for  air  temperatures  between  -173  and  277°C. 

Friedlander  (1977,  p.  105)  gave  a  modified  form  of  the  drag  coefficient,  eq  77,  that  is 
more  accurate  at  higher  Reynolds  numbers 


CD  =^a-(l  +  0.158  Re  2/3 ) 
u  r 


(79) 


Here  the  Reynolds  number  is  defined  as 


Re  =ZJJI 
v. 


(80) 
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Substituting  eq  79  into  eq  76  yields  an  equation  in  u( 


9v,[l  +  0.158(2  r  uf/va)2/3]  'P*  '  (81) 

It  is  routine  to  solve  this  equation  iteratively  for  uf  as  a  function  of  r0  by  Newton's  meth¬ 
od. 

Figure  15  shows  typical  values  of  uf.  Notice,  modifying  the  drag  coefficient  according 
to  eq  79  is  important  for  the  largest  droplets;  the  curvature  in  the  figure  at  the  largest  ra¬ 
dii  indicates  that  uf  is  deviating  from  the  proportionality  to  r2  represented  at  the  smaller 
radii. 

In  still  air,  jet  droplets  generally  are  projected  to  no  more  than  20  cm  above  the  water 
surface  (Blanchard  1963,  Bortkovskii  1987).  The  turbulent  air  over  the  ocean  will  surely 
lift  some  of  these  to  greater  heights;  it  will  undoubtedly  suspend  most  of  the  film  drop¬ 
lets  and  may  even  carry  the  large  spume  droplets  upward.  The  spray  droplet  profiles  that 
de  Leeuw  (1986, 1987)  measured  over  the  open  ocean  are  evidence  of  this  turbulent  mix¬ 
ing.  His  profiles  show  very  little  decrease  in  droplet  concentration  from  the  surface  to 
heights  of  10-20  m  for  droplets  with  radii  less  than  50  pm,  the  maximum  size  he  could 
observe.  As  it  turns  out,  the  fate  of  spray  droplets  is  one  of  the  unanswered  fundamental 
questions  in  this  field.  What  is  the  percentage  of  droplets  that  immediately  fall  back  to 
the  surface,  and  what  is  the  percentage  of  droplets  that  remain  suspended  by  the  turbu- 


Figure  15.  Values  of  the  terminal  fall  speed  computed 
from  eq  81 .  Environmental  conditions  are  Tfl  =  -10X1, 
S  =  34  %c  and  P  =  1000  hPa. 
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lence  and,  thus,  approach  their  full  thermodynamic  potential?  Turbulence  modeling  of 
the  marine  surface  layer  like  that  begun  by  Fairall  et  al.  (1987)  is  vital  to  addressing  these 
questions. 

It  is  therefore  an  oversimplification  to  base  spray  droplet  residence  times  only  on  the 
terminal  fall  speed.  If  the  fall  speed  is  less  than  typical  turbulent  vertical  velocity  fluctua¬ 
tions — nominally  0.2-0.3  m/s  in  a  10-m/s  wind — the  turbulence  may  suspend  the  parti¬ 
cles  indefinitely.  According  to  Figure  15,  droplets  with  radii  less  than  about  50  p m  have 
terminal  fall  speeds  less  than  this  limit.  Only  for  the  largest  droplets — for  which  the  fall 
speed  is  much  greater  than  the  turbulent  velocity  fluctuations — can  we  be  certain  that  the 
physics  contained  in  eq  75  is  sufficient. 

The  terminal  fall  speed,  nevertheless,  does  give  us  an  intuitive  feel  for  the  atmospheric 
residence  times.  In  Figure  16  I  thus  plot  a  time  constant  or  residence  time  rf,  which  is  the 
time  required  for  a  droplet  of  radius  r0  to  fall  1  m  in  still  air.  That  is. 


(82) 


For  comparison,  in  the  figure  I  also  plot  the  tT  values  and  two  sets  of  the  tr  values 
from  Figure  13b.  We  see  that  many  of  the  spray  droplets — namely,  those  with  radii  great¬ 
er  than  about  20  jum — will  probably  fall  back  into  the  sea  before  reaching  their  equilibri¬ 
um  radius.  Because  of  the  rapidity  of  the  thermal  exchange,  however,  all  but  the  very 
largest  spray  droplets  should  reach  thermal  equilibrium  before  falling  back  to  the  sur¬ 
face. 


Figure  16.  Time  constants  for  thermal  (tt)  and  mois¬ 
ture  (xr)  exchange  and  the  time  required  for  a  droplet 
to  fall  1  m  in  still  air  (if).  Environmental  conditions 
are  T W  =  0°C,  T(I  =  -10°C,  S  =  34  %c  and  P  =  1000 
hPa. 
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CONCLUSIONS 


The  equations  that  I  have  presented  constitute  a  powerful  model  for  treating  many 
problems  related  to  the  evolution  of  sea  spray  droplets.  I  have  shown  a  few  applications 
here,  such  as  modeling  how  droplets  approach  thermal  and  moisture  equilibrium  and 
finding  time  constants  characterizing  these  processes.  An  important  conclusion  of  these 
modeling  runs  is  the  disparity  between  the  time  constants  for  thermal  (t,-)  and  moisture 
(tr)  equilibrium.  Although  both  Zj  and  rr  increase  with  the  initial  spray  droplet  radius 
roughly  as  r02,  at  a  given  radius  xr  is  always  at  least  three  orders  of  magnitude  larger  than 
tt  for  humidities  typical  over  the  high  latitude  ocean.  Consequently,  the  ambient  relative 
humidity  has  negligible  effect  on  the  thermal  evolution  of  spray  droplets;  and  the  air-sea 
temperature  difference  has  minimal  impact  on  their  moisture  evolution. 

Although  all  the  computations  I  discussed  were  for  a  steady-state  environment,  the 
model  can  easily  handle  droplet  evolution  in  a  changing  environment.  This  capability 
will  allow  us,  in  the  future,  to  treat  not  just  single,  isolated  spray  droplets  but  also  a 
cloud  of  interacting  droplets.  The  groundwork  is  now  in  place  for  our  research  on  the 
role  of  sea  spray  in  the  air-sea  exchange  of  heat  and  moisture. 
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